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IKE  PRODUCTION  OF  TURBULENCE* 

By  W.  Tollmien 

It  is  known  that  for  geometrically  similar  boundary  surfaces 
the  form  of  flow  of  liquids  and  gases,  in  so  far  as  the  effect 
of  compressibility  can  be  disregarded,  is  determined  only  by  the 
Reynolds  Number,  that  nondimensi  onal  or  absolute  quantity  which 
is  composed  of  a  characteristic  velocity  and  distance  and  the 
kinematic  viscosity.  The  most  conspicuous  and  iv^-ortant  change 
in  the  form  of  flow  with  changing  Reynolds  Number  is  the  sharp: 
transition  from  the  laminar  to  the:  turbulent  state  of  flow.  The 
Reynolds  dumber  for  this  transition  has  been  determined  from 
numerous  observations  in  many  cases;  for  example,  in  the  case 
of  the  flow  through  tubes  or  about  a  sphere.  Still  the  summary 
announcement  of  this  "critical*  Reynolds  Number  is  physically 
incomplete  since  it  depends  largely  on  the  disturbances  of  the 
original  laminar  flow.  From  the  disturbances  we  now  have ,  at 
the  worst,  the  statement  of  their  point  or  origin,  e»g« ,  in  the 
case  of  the  flow  through,  tubes,  whether  the  disturbances  proceed, 
from  the  inflow  or  farther  downstream  from  the  wells,  where  they 
apparently  have  less  effeer  than  in  the  former  case.  Beyond 
this  point,  however,  the  disturbances  hate  not  ret  been  exaeri- 

r.entelly  defined  and  classified.  'Hence  the  experiments  leave  it 

*“Ueber  die  Ente.tehung  der  Tuxbulenx. H  From  I-achrichxen  dor 
Oeselleohaft  der  Vi s  eenschaf  ter.  z\k  Gottingen  (Report  I  Of  the 
G^t tinker  Scientific  Society),  the.rati sch-Pky si ksl i sche  Elacse, 

1S?9,  pp.  ?l~4f . 
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uncertain  in  what  sense  one  can  speak  accurately  of  a  oritical 

Reynolds  lumber  as  an  upper  limit  for  the  laminar  flow,  and 

> 

whether  this  flow  can  be  maintained  up  to  a  large  Reynolds  Hum¬ 
ber,  by  continually  reducing  the  magnitude  of  the  disturbances. 

n 

We  might  hope  to  obtain  an.  accurate  insight  into  the  pro¬ 
duction  of  turbulence  by  a  photographic  investigation  of  the 
changes  in  the  flow*  Reynolds1  method  has  been  much  used,  ac¬ 
cording  to  which  the  flow  in  glass  tubes  is  observed  with  the  at 
of  a  colored  thread  suspended  in  the  middle  of  the  tube.  Prandt 
sought  to  obtain  a  clearer  idea  of  the  process  in  the  smooth 
flow  of  Water  in  an  open  channel  by  photographing  aluminum  pow- 
der  strewn  on  the  surface.*  By  this  means  he  succeeded  in  de¬ 
termining  approximately  the  location  of  the  flow  transformation, 
e.g. ,  in  the  flow  through  a.  tube,  mostly  in  the  so-called  start- 


to  those  obtained  by  pressure  and  velocity  measurements 


only  in  rough  outlines.  Hitherto  the  pro  auction  of 


turbulence  could  not 


In  any  case,  according  tc  the  customary  view,  the  produc 
of  turbulence  begins  with  the  laminar  flow  becoming  unst 


obtain  a  theoretical  conception  -of 
a  level  laminar  velocity  distribution 


bate  perpend! cu- 

:  dm),  p.  4?: 
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lar  to  the  direction  of  flow)  and  then  to  subject  this  Telocity'* 
profile  to  two-dimensional  disturbances  *  She  first  profile  to 
he  treated  in  this  way  by  the  method  of  small  oscillations  was 
that  of  the  “Couette  flow”  with  linear  dependence  on  the  veloc¬ 
ity  of  the  transverse  ooordinate.  The  result,  obtained  by  Kopf* 
and  liisee**,  gave  no  indication  of  the  production  of  turbulence, 
but  a  disappearance  of  all  disturbances,  and  consequent  perma¬ 
nent  stability  of  the  assumed  laminar  flow*  The  simplification 
of  the  disturbances  was  first  obtained  by  Prandtl  (See  footnote* 
page  3)  and  Tietjens***  for  certain  profiles.  These  investiga¬ 
tors,  however,  obtained  the  physically  surprising  result  that, 
within  the  range  of  the  Reynolds  Numbers  for  which  the  approxi- 

s: 


i 

I 


while  it  was  to  be  expected  that ,  belovr  a  certain  Reynolds  Num¬ 
ber,  the  laminar  flow  would  be  stable  despite  all  disturbances. 

In  the  Prandtl-Tietjens  investigations,  for  mathematical  sim¬ 
plicity  according  to  a  suggestion  of  Rayleigh,  profiles  were 
assumed  for  the  fundamental  flow,  which  were  composed  of  straight 
sections.  The  universal  disappearance  of  the  curvature  of  the 
profile  of  the  fundamental  flow  indicates,  however,  in  this  Case, 
such  an  essential  physical  assumption  that,  by  allowing  for  the 
curvature,  a  modification  of  the  Prandtl- Tie  tie  in  result  might 


be  expected  at  the  ; 

*L. 

Her f ,  Ann.  d. 
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Weber 

Festschrift, 

***c 

Piet tens,  2. 

300. 

m 


4 


K.A.C.A.  Technical  Memorandum  Ho.  609 


what  follows*  therefore,  profiles  axe  used  which  do  not  have  a 
universally  disappearing  curvature. 


This  introduction  is,  of  course,  not  in  the  least  intended 


give  a 


survey  of  the  extensive  liierature*on 


the  production  of  turbulence.  It  is  only  intended  to  indicate 
the  connection  in  which  the  present  work  was  undertaken.  I  am 
greatly  indebted  to  Dr.  Prandtl  for  suggesting  this  investiga¬ 
tion  and  for  his  encouragement  and  interest. 


1.  General  Disturbance  Equation 


The  principal  flow 


,  y  being  the 


coordinate  perpendicular  to  the  direction  of  flow  in  the  t 
axis.  The  disturbances  are  resolved  according  to  Fourier,  so 
that  the  flow  function  of  a  given  Fourier  term  (y)  e* (h-x-PO  » 
<Ry)  .  Then  tt  represents  the  spatial  cycle,  the  real 

part  of  B  represents  the  periodic  oycle,  and  the  real  part  of 
c  represents  the  phasic  velocity  of  the  giver,  disturbance; 

The  imaginary  part  of  P  indicates  the  increasing  or  damping  of 
the  disturbance  according  to  whether  it  is  positive  or  negative. 
According  to  the  method  of  small  vibrations,  only  the  linear 
terms  in  V(y)  are  always  considered.  Then,  by  introduction 
into  the  favier- Stokes  differential  ecu? ti ons  for  viscous  fluids, 
eliuihatior.  of  the  press-ore  and  general  introduction  of  absolute 
quantities,  we  obtain  the  disturbance  equation 


(V  -  c  *  -  d*  c) 


W  --  v  -4  V  ;  ;  V'  V;tT  ■»  L-  <x  tT 


Ov  4  C*  4'f 


R.A.C.A.  Technical  Hemoxabdum  Ko.  609  5 

in  which.  R  is  the  Reynolds  Humber .  The  reference  length,  with 
which  a  characteristic  velocity  of  the  profile  (e.g.,  the  maximum 
Xjm)  multiplied  after  division  hy  the  kinematic  viscosity  yields 
v  R,  is  so  selected  as  a  dimension  of  the  profile,  that  the  de¬ 
ductions  of  the  absolute  U  in  equation  (!)  are  of  the  order 
of  magnitude  of  1.  The  disappearance  of  both  components  (<p  and 
q>*)  of  the  disturbance  is  a  marginal  condition  at  the  boundary 
walls,  i . e. ,  the  limits  of  TJ(y)  by  definition. 

At  first  tee  will  assume  f  and  c  to  be  real,  and  will 
seek  the  points  of  transition  from  stability  to  instability. 

For  this  problem  two  statements  can  now  be  .made  concerning  c 
and  a*  If  we  write  the  disturbance  equation  symbolically 
t>(<p)  =  0  and  integrate  the  expressions  f  L(<P)  -  <P  L($)  and 
?  L(<P )  +  cp  L(  (the  conjugated  complex  quantities  being  over- 
lined)  between  the  boundaries  with  allowance  for  the  boundary 
conditions,  we  then  obtain,  in  the  Simple  and  well-known  manner,' 
c  <  for  Un  §  0  and  also  a  <<  &-**  Just  the  circumstance 

that  we  must  accordingly  assume  a  zero  position  of  U  -  c,  ne¬ 
cessitates  special  mathematical  expressions  and  is  also  accotjva- 
nied  by  a  physically  interesting  behavior  of  the  disturbances. 


•footnote  from  page  4:  * 

T .  Foe  then,  Z.  f.  anrew.  Hath.  u.  Keen,  1,  page  135,  1921; 

L.  Schiller,  Rhys.  Z~  26,  page  566,  1935.  \ 

*  *  Sole  erg,  "The  Turbulence  Problem,1’  Proc.  of  the  First  Interna¬ 
tional  Congress  for  Applied  Mechanics,  Delft,  1924,  page  38$. 

Roether ,  ”  Zur  asympccti sober  B eh and lung  der  stationSf en  Lo sun¬ 
der  i;  Curbulenzororlem, R  Z.  f.  ansew.  -Hath.  u.  Hech. ,  Vol.  6, 
1926,  cafe  232. 
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II,  Frictionless  Solutions  of  the  Disturbance  Equation 


In  general  a  R  is  as  sinned  to  be  very  large.  It 
to  endesCvor  to  represent  the  solutions  of  the 


is  natur- 


equaticn  in  such  manner  that  the  terms  multiplied  by 

1 


^  can  be  omitted,  so  that  they  axe  equivalent 
to  solutions  of  the  frictionless  disturbance  equation # 

(U  -  c)  («P«  ~  as  q>)  -  U«  <P  =  0 

We  will  consider  the  solutions  of  this  frictionless  disturbance 
equation  without  at  first  going  further  into  their  connection 
with  the  solutions  of  the  general  disturbance  equation  (for 
large  aR). 

At  the  position  U  -  c,  the  differential  equation 

U* 


point* 


a 


we  can  easily  establish,  ac 


U  ~  c 

in  this  position,  provided  Un  does  hot  vanish  at  just  this 
point,  we  have  a  nonessential  singularity  (point  of  precision) 
under  consideration,  a 

regarding  linear  um twin, +?>.*.  cuuavi,w«D> 

■xng  series  of  integrals  from  the  point  V  -  «  =  0.  At 
this  point  let  y  »  0  and  y  be  positive  for  N  -  c  >  0.  Ac¬ 
cordingly  the  values  at  this  point  are  indlca  tedv.by  the  Subindex 
C.  mere  are  then  two  linear  independent  solutions  of  equation 
(2),  whet  the  potential  series  in  y  are  designated  by  | %'W) 
and  Pg ( y ) ,  and  when  the  constant  determined  by  the  differen¬ 
tial  equation  is  designated  by  A. 


fb  - ?  ******* 
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%.  =  y  Pj 


y  +  . . . 


.  V  u»‘ 

%  ~  %  (y)  +  A  %  log  y  -  1  +  ...  +  ~_£  y  log  y  +  ... 

U*o 

.B  establishes  the  fact  once  for  ail  that  log  y  is  real 


Before  drawing  conclusions  from  the  analytical  character 
of  <Pj  and  %,  these  should  be  determined  for  a  special  case 
of  U  and  indeed  for  a  parabolic  velocity  distribution.  Ac¬ 
cordingly,  U  -  c  «  y  {2a  -  y),  so  that  U*  =  0  for  y  «  a. 
Then  U  -  1,  if  we  put  a  *  J  1  -  c,  and  17  =  0  for  y  - 


a  -  1.  We  expediently  introduce  ya  «  y/a  and  a*  «  a  a,  so 
that  ys  -  1  at  the  vertex  of  the  parabola,  and  we .  then  obtain 

yx  (2  -  y, )  <0«  -  a if?  q>)  +  2  <p  =  0 


as  the  differential  equation  in  the  new  variables.  We  then 


obtain  the  series  developments 


?*■  “  +  ^  y  *  +  yf  +  , .  , : 


whereby 


n  (r  -  3}  ar,_j  + 


ai£  an-3 


n  In 


ai  “  I*  ;%  ~  =  •%—:>  a*  -  - 


w hereby 


■ffe  s  1  t  tj  y.  t  \  yxs  +  bs  y;s  ■+  . . .  -  log  y 

V  ~  *'n~5)br,_.yf3ab  Zr.-~.i-~S  W342{2n-l)ar)-(2n-3)aT^. 


5»  ■  —  7b 
£  i.  V 


'5*  bh  it  ^  -i*  ..  -,Z  —  ^  7-“  ££  -*4  ~Yn 
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In-  the  general  consideration;  of  <Pa  and  '  fs  we  first 
find  that  fv  remains  regular  even  at  the  singular  ;point*.:  while 
fs  has  a  logarithmically  infinite  derivation  at  that  point. 

If  we  remain,  at  least  in  so  far  as  c??  is  conoefned,;:  a;t;;ioae::;:;:;:; 
distance  from  the  singular  point  U  m  c,  then  :Cpj  .  and:  9>-:  ;ap-x; 
proximate  the  complete  disturbance  equation  (l)  for  a  very  large 
O  R  :and  all  the  closer,  ■: the  greater  a  H  le.*- : -.er rorj: ln£ ••: 

and  $£  is  small  and  of  the  order  (c.  B)  \  On  account  of 
the  infinity  of  the  derivatives  of  %,  this  conclusion  |;i  ob¬ 
viously  no  longer  justified  $2i’ ' Hi : *£3LH||g I ?o||! [ ^ i- 
critiosl  point.  We  must  therefore : discuss ; : the :  general  disturb- 
: ance  ■  equation  (l)  in  the  vicinity ;  of  *  V  * : c.  This  point  is 
^physi c?  lly  character ized  by : :  the  fac t :■  that : due .V; 

Of  the :  main :  veloci ty  and: :  the  :  phasic ■ : velocity,  any  fluid  par t  i~  : 
cie  always'  remains'  in  the  seme  pressure '  field,  if  the  slight:',  -a 


,e  frictiorless  solutions  in  a  more  formal 

:^8iM^::::::dhtdgrali  ]'\M  sturfe; :  •  ...• ; 

aside  from;  tker-cxii leal,  point  V  *»  C:»  re~; 
ir  derivatives  for  any  increasing:-,  vaiue.e-:- 


■t<  -i 


;  >H- 
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tions  of  the  Disturbance  Equation 
in  the  Vicinity  of  U  =  o 


uiiC  investigation  oi  tne  solutions  in  the  vicinity  of 
the  critical  point  U  -  c,  we  take  a  small  region  about  this 
point,  the  “transition  strip,”  in  which  V  -  c  can  be  replaced 
with  sufficient  accuracy  by  U*D  y  and  likewise  U“  by  u«0. 
te  nos  introduce  y  =  (a  R  C*0)-l/3r,  ,  £  t), '  whereby  r  my  ims 
very  great  values,,  due  to  the  smallness  of  c ,  oven  at  small 

values  of  y.  From  equation  (l)  we  then  obtain  the  differential 
equation  in  n  : 


vre  obtain 


n,  Solutions  of  the  dis- 
small  c)  in  the  vicinity 


equation  for  large  aR  (i,e 


equation 


point 


proximate  solutions 


follows 


critical  noint 


on  ration  o:  v£  c  ,  introduce  into  the 
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U«  .  .  .V  . 

term  e  <P  and  calculate  from  differential  equation  (4)  a 
correction  which  we  will  designate  by  c  vsl.  The  differential 
equation  fox  fsl  Is  then 

-  i  0>™  =  T,  <p.„  _  g^S  . 

Since  it  may  be  assumed  that  the  friotionless  solution  is  appli- 

oable  at  some  distance  from  the  critical  point,  we  will  fry  to 

obtain  a  solution  of  equation  (5)  which  is  based  outwardly  on 

the  corresponding  frictionless  solution*.  In  the  vicinity  of  the 

Critical  point  (i.e. ,  for  small  Jyj  )  a  friotionless  calcula- 

tion  yields  <p|x  =  rr~  _  .  This  value  mast  serve  therefore  as 

u  o  1 

the  basis  for  the  solution  Of  equation  (5)  outwardly;  for  exam¬ 
ple,  for  large  {rj{  * 

1%  by  way  of  abbreviation.  We  designate  by  kx/3  and 

J^S)  . 

:Ki/i  the  Kanhel  functions  of  the  first  end  second  form  with 
the  eubir-dex  1/3  and  the  argument  |-  (in)3 /a  :  £  *" 

then  the  solutions  of  the  homogeneous  differential  equation 
( 5) ,  namely 

S  p}  '  an^ 

r*l/3  a-“‘~  *  “1/3  VP/ 

yield  a  solution  of  the  inhomogeneous  differential  equation  (5} 

H  r,l/  fT  dt  t'  r  1  »r\  fr)\ 

u*A  6  in/3  J  •“i/s  dr  -  r*x/s  J  ri  r.i/3dn  (7)  \ 


.  h\ 

h;;5  dt  - 

v(l)  rf>  rl/3T.(*)  /„* 

hi/s  J  r<  )  (?)  | 

tailor  of  the 

H-mkel  functions  p  We 

in  the  lower 

h  half-r lane  (which 

laxity  of  one 

solution'-,-  s-:  that  . 
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0  *  axg  T)  ~  -  tt.  Moreover,  ^  «  arg  C  (in)  3  £  -  For  the 
connection  with  the  fricticnless  solution,  we  utilize  the  be¬ 
havior  of  equation  (?)  fox  large  M  .  At  the  junction  points 

The  lower 
a  point 

semicircle  in  the  lower  half -plane.  We  replace  the 
functions  by  their  asymptotic  formulas  for  large  values 


h  =  *  r.  Gn  the  circle  with  this  radius 
limit  of  both  integrals  of  equation  (?)  < 


notes  the  exponential  function  of  the  bracketed  expression. 

For  the  following,  we  need  equation  (7)  on  the  semicircle  in  the 
lower  Tj  half-plane  with  r.  The  first  integral  in  equation  (?) 
there  attains,  in  the  first  asymptotic  approximation,  the  value 


Memorandum  No 


Thus  we  obtain  the  desired  <p« 


equation 


—  Tai  ^  iiop  Y'S1  according  to  %  whereby, 

ae  the  integration  method,  we  again  use  the  semicircle  on  which 
is  now  known.  This  is 


in  a  similar'  manner  it  is  finally  found  that  the  frictionles 

solution  on  the  positive  edge  (+r)  of  the  transition  strip 
becomes 


known  to  be 


uxior.  in  5 


v':-' 
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at  the  critical  point ,  which  consisted  in  the  introduction  of  an 
imaginary  portion  and  therefore  of  a  phasic  jump  in  the  *  com¬ 
ponent  of  the  disturbance.  The  above  calculations  show  that  the 
determination  of  the  transitional  substitution  is  equivalent  to 
the  determination  of  the  physically  real  branch  of  the  logarithm 
in  9  s' 

statement  may  be  made  regarding  the  width  of 
strip,  for  which  we  deduced  the  <PEi  solution, 
r  must  be  so  large  that  the  asymptotic  approximations  for  large 
r,  as  used  in  the  above  calculation,  are  applicable  up  to  an 
error  which  is  admitted  on  account  of  the  only  approximate  appli¬ 
cability  of  the  frictionless  solutions  outside  of  the  transition¬ 
al  strip.* 

Calculation  from  the  thus-determined  r  yields  2  e  p  == 

2  (a  H  U*p)  r,  as  the  width  of  the  transition  strip  in  the 
original  y  coordinate.  Hence,  as  soon  as  aR  is  large  enough, 
the  frictionless  cp”  x,  due  to  the  narrowness  of  the  transition 
Strip  which  approaches  0  at  the  limit,  actually  acquires ,  on 
its  boundaries,  the  values  We  have  above  put  ir.  the  boundary  con¬ 
ditions  .  Likewise  ?2,  in  the  above  integra tier,  method,  suffi¬ 
ciently  approximates  the  value  1  Kh.ch  we  assumed  for  ,  thus 
verifying  the  asymptotic  applicability  of  our  formulas  for  large 
Q  R,  Moreover,  in  Section  VI.  3,  reference  is  made  to  a  somewhat 
-rent  and,  expander  derivative  cf  the  transitional  substitution 


the  asymptotic  development  cf  fos  large  the  term 

-e  —  U,!r  I  y  if  c  i  *  1 

following  d  iE  pr ,  fbcm  which  the  indicated  error  car 

be  estimate*.0  1  w  0  r' 


:+  i- 


mmm 

wmM 
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Figure  1  gives  for  UM0/U»o  =  1  a  graphically  obtained 
qualitative  picture  of  the  behavior  of  V  21  in  the  transition 
strip,  divided  into  a  real  part  R(q>  *S1 )  and  an  imaginary  part 


The  differential  equations  (4)  and  (5)  give  us  directly 
the  second  pair  of  solutions,  which  follow  from  the  general  dis¬ 
turbance  equation  for  large  aR.  Along  with  the  solution  <jp» 

«  X 

of  equation  (5),  which  we  have  just  considered,  there  are  obvi-. 
ously  possible  also  the  solutions  (equation  6)  of  the  correspond 
ing  homogeneous  differential  eauation. 


written 


follows 


The  error  in  <P3  is  of  the  order  of  magnitude  of 
According  to  their  rate  of  variation  cpa  and  ?£ 
designated  as  *6lor"  and  <?-  and  as  ’’rapid.1* 


iiis 
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■'  i*  ’ 

IV.  Solutions  of  the  Disturbance  Equation 
on  the  fall  for  c  »  v 


the  distance  from  the  wall  is 


the  solu¬ 
tions  of  this  differential  equation  have  thelforin  %  • 

}  A  **  w 

'W  ■*.  -  ^  ..  ^  r  £  ^  .  W-  v* 

te  thus  obtain  both  pairs  of  solutions  in  an  especially  clear' 
fanner.  We  have,  namely,  *  61  *  >"1 

*\a  !  /  T  '  V  *  , — , 

k  =4j-  iaRc  +  3a8  ±  J  -  a*  R2  c*  ~  4  i  a  ft  > 


One  k  accordingly  approximates  a2 - -  4  i  ,— -y,„,  and  gives 

45  a  Rc3 

us,  as  a  result  of  the  exclusive  consideration  of  the  largest 
terms  with  large  OR,  the  pair  of  frictlonlees  solutions,  i.e* 
due  to  the  dependence  of  the  coefficients  of  the  differential 
equation  of  yw,  only  the  beginning  of  their  development  from 
the  wall  out.  approximately 


tf.A.C.A.  Technical  Memorandum  flo.  609 


The  solution  Corresponding  to  <p: 
a  certain  constant  factor 


{-l+^/aRcy^ 


i.e.,  it  diminishes  very  rapidly,  sp  that  the  variation  o?  the.h! 
coefficients  of  the  differential  equation  does  not  enter. Into 
consideration*  It  should  he  here  noted  that  equation  {10}  can^ 

,  -  'i'jr  JL, 

also  he  derived  from  the  general  equation  (8)  for  <Pa  by 

cialization  for  negatively  large  tj. 

Tiet^ens  derived  <p3  from  a  separate  boundary-layer  ^qua-H 

-  il- 

tion,  while  we  Simply  derived  all  solutions  from  the  general 
turbance  equation.  <p3  can  he  conceived  as  a  sort  of  boundary*- J 
layer  solution  on  the  wall  or,  in  other  words,  of  outer  frictleii 
al-layer  solution.  In  the  case  U*  4  0,  however,  at  the 

;  .  J 

ical  point  where,  according  to  the  fxictionleBS  calculation,  1; 
one  of  the  disturbance  components  would  become  infinite,  a  see-  . 

. .  ,  ..  .  :  m 

oiid  inner  fri'ctional-layer  solution  occurred  which  was  derived  % 
as  <f>zx  in  Section  III.  In  the  above-mentioned  special  case  of 
greater  c,  the  two  frictional  layers  also  separate  in  spao'e**^ 


in  our  investigation,  as  Compared  with  that  of  Tietjens,  r 
in  the  appearance  of  the  inner  frictional  layer. 


*v 
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V*  Formulation  of  the  Boundary  Conditions 


Having  obtained  an  insight  into  the  character  of  the  four 


integrals  of  the  disturbance  equation,  we  will  now  revert  to 
boundary-value  problem  and  formulate  the  boundary  conditions 


,  which  increases  from  the 


We  first  take  a  profile 


value  U  »  0  on  the  wall  to  a  maximum  value,  which  then  remains 
constant  for  any  desired  time.  This,  for  example,  holds  very 
true  for  the  important  profile  of  the  flow  along  a  plate.  Since 
for  positive  y  increases  beyond  all  limits,  it  does  not 
appear  in  the  solution.  For  TJ  =  constant,  the  flibtioniess 


solutions  are  e~ 


is,  we  can  use 


due  to  the  necessary  limitation. *  Hence  cp>/<p  =  ~  a  in  this 
field,  whereby  we  have  accordingly  obtained  a  simple  boundary 
condition  for  the  initial  point  of  the  field  with  constant  17, 
which  will  be  indicated  in  future  by  the  subindex  m.  Since  cp3 
has  already  disappeared,  the  condition  is  expressed  by  the  equa- 


4  a  ^im'  4  1  r  m  4  c  :4?  0 

or,  in  an  abbreviated  form,  by 

C,  %m  4  %  %  -  C 


(lit) 


0  an;  <p 1  must  how  vanish  at  the  wall.  If  we  de  si  gnate 
the  v&iuet  there  by  the  subindex  •  w  and  simultaneously  apply 


e4a>*  and  cannot  balance  each  other,  because  c  <  Jc.  K 

feocordinr  to  a  remark  in  Section  1 . 
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substitution  in  %w  and  in  <p»aw,  we  obtain 

••  •  ■  ’  '■  •• :  •  ■  .  •:  •  :  . .  ’  •  ;  /  :■■■■' ■.■■v 

$1  ^iw  *  Cj  ^sw  +  C3  <P3W  »  0  (1B&) 


and 

or 

The  ratio 


Pi  +  0* 


<P*aw  +  cs 


n 


(i3by 


■<p 


Cl  %,  w  +  9».w  -  (0i  <%*  +  <PV*>  -  0 

3|f 


has  already  been  calculated  in  the  work  of  Tint- 

T  aw 


jens  and  can  be  expressed  by  the  symbol  B.  It  then  becomes 

> 

^3  XT 

sir~  ^  w  c  ^ 
y  aw 

D  is  here  introduced  as  a  function  of  z0,  which,  in  the  % 

•  .  ’  :'"v  f 

system  of  coordinates  (Section  III),  is  the  measured  distance 
of  the  wall  from  the  critical  point,  which  can  be  designated  by 
hence 


o  1 

U*r>  € 


The  following  table  gives  the  values  according  to  Tietjens 


TV 

*W 

0 

D(V 

0.702  -  0.425  1 

Uw 

-3.0 

D(hw) 

1.400  +  0.515  i 

— 0. 5 

0.785  -  0.411  t 

-3.5 

1.180  * 

1*130  1 

-1*0 

C.  920  -  0.389  i 

-’4.0 

0.460  + 

1.250  1 

-1.5 

1.043  -  0.297  i 

-4.5 

-0.0405  + 

C .  6080  i 

-2.0 

1.206  -  0.147  i 

-5.0 

0.0057  + 

0.3645  i 

-2.5 

1 . 357  +  0. 106  i 

-5.5 

0.1913  + 

0.2393  i 

The  two 

conditions  (lla)  and 

(13c)  v 

ield  a  det 

eminent  relation 

which  car.  be  calculated  from 

Svr  /T  ’  ; 

*  We  obt 

sin  thereby 
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&  m  _  &  m 

2vj  v2m  Tiw 


''sia  Tiw 


t  = 


0>o  ^tra  ~  ^iw 

c  ^irn^sw  -  ^a»  ^vw 


fhe  right  Bide  of  equation' (15)  is  known,  when  we  write  the 
parameter b  c  and  a  before  it.  It  is  quite  generally  desig¬ 
nated  by  E ( c :  c) .  The  left  side  is  expressed  as  a  function  of 
r^m  We  now  construct  a  polar  curve  from - with  its  real 

part  as  abscissa  and  its  imaginary  part  as  ordinate  (Fig.  3). 

—■  •  * — 

We  then  select  a  definite  c  and  plot  E  as  a  function  of  <x 


in  the  polar  curve,  whereupon  we  determine  the  intersection 

D  (rjw) 

points  with  the  -  polar  diagram  and  the  %  correspond¬ 

ing  to  the  n transni ssi on  point11  (Section  I) From  %  we  can 
then  determine  c  according  to  equation  (14)  and  consequently 


the  last  parameter  B. 

Regarding  D  it  should  be  noted  " 

eW* 

large  rt  %  it  is  convertea  into-  p";;*: 
pushed  either  through  an  asymptotic  c< 
pres  si  or.  for  as  given  ir.  equation 


y  for  a 


—  .  This  can  be  accom- 


:  ex~ 


derivea  equation  (10)  ir. 


expressed  by  riW. 


or  from 


c  is  t 


For  a  symmetrical  !■  profile  we  obtain  a  simple  form  of 
the  boundary  condition  from  the  fact  that  in  this  case  s  Cain 
be  resolved  into  an.  exren  and  an  uneven  part,  each  of  which  must 


•  -  HHiWljT':-* 


Memorandum  No 


of  itself  satisfy  the  disturbance  equation  and  the  boundary  Con¬ 
ditions.  Hence,  we  have  in.  the  middle  j  since  <p*  cannot  be  in¬ 
troduced  here  and  q>3  has  already  vanished,  the  condition 


for  the  uneven  V,  which  takes  the  place  of  equation  (ll), 
while  the  condition  at  the  wall  (equation  12)  is  maintained. 

Tietjens  expressed  the  one  boundary  condition  on  the  wall, 
and  indeed  for  the  y  component  according  to  a  suggestion  by 
Prahdtl  in  a  different  form  from  ours.  Since  doubt  has  occa¬ 
sionally  been  expressed  regarding  the  correctness  of  this  bound' 

ary  condition,  we  hereby  confirm  its  agreement  with  our  own  con- 

•  _ 

•j  \ 

elusions. We  will  describe  the  Tietjens  method  formally  and 
refer  to  the  work  itself  for  its  physical  justification,  Tiet- 
jens  doer  not  let  his  frictional-layer  solution  for  H  *=  +  •, 
which  would  correspond  to  our  %  and  be  designated  by  qy* 

{ proper tional  tc  Tietjens  *  v*) ,  go  toward  0  like  ours,  but  to 
ward  &  constant  value  *¥s  V  (proportional  to  v*),  so  that 

~  $3  *  %<?y  for  <?**'*•,  like  ?-*3  ,  is  also  put  equal  to  : 
Then  the  constant  is  so  determined  by  Tietjens  that 


one  for  the  wall,  is 
ed  by  numbers. 
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to  v  that  is. 


Cx  <Ptw  * 


■'While  we  above  put 


Ci  <PiW  + 


+  O3  <P3W  &  0 


Both  conditions  are 
tion  between  %£ 


therefore  identical  according  to  the  rela- 


VI.  Equilibrium  of 


1.  The  above  method  of  calculating  will  be  applied  first 
to  the  velocity  profile  which  develops  at  large  Reynolds  lum¬ 
bers  on  a  plate  immerged  parallel  to  a  flow  of  infinite  extent, 
the  leading  edge  of  the  plate  being  perpendicular  to  the  direc¬ 


tion  of  flow.  We  choose  this  particular  velocity  distribution, 
because  experiments  on  its  stability  are  in  progress  and  because 
the  assumptions  for  the  accuracy  of  our  approximation  are  here 


especially  favorable.  This  profile  was  calculated  by  Blasius 
according  to  Pr and tl *s  theorems^  and  is  Characterized  especial¬ 
ly  by  the  fact  that  it  begins  at  the  wall  with  a  vanishing  cur¬ 
vature  and  joins  the  undisturbed  velocity  outwardly  with  a  sharp 


to  the  asymptotic  be! 


\e  error 


integral  (Fig.  2). 


our  purposes 


usef 


representation  of  the  profile.  According  tc  the  above*  the  sim- 
- - —••••••••• — — £ - - - - - - - - - -  - - — — ■  ■  . . 

;L.  Frandtl,  ''Ueber  Fldssigkei tebewegung  bed  eehr  kleiner  R:ei- 

cung,'»  Vex'd,  d.  Ill  Internet.  Hath .  Kongresses ,  Felder:  erg,  2924; 

n | printed  in  Vier  Abhend  lungers  znr  Hy drody nafeib  und  Aerciynsmik, 

lottingen,  1927.  For  translation,  sec  F . A . C . A .  Technical  Heme- 

rap^-nrr  Va  AC,? 
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ples±  way  is  to  let  the  velocity  distribution  begin  at  the  wall 
with  a  straight  section  and  to  replace  the  asymptotic  junction 
with  the  undisturbed  velocity  by  a  finite  junction  (at  M  in  Pig. 
2),  i.e.,  to  continue  the  straight  section  by  a  parabolic  sec¬ 
tion.  If  we  designate  by  ym  the  distance  from  the  vertex  of 
the  parabola,  we  then  allow,  toward  the  one  side,  U  constant 
to  continue  with  its  maxi  cram  value  (l).  Toward  the  other  side 
U  =  1  -  ya£ ,  if  we  correspondingly  select  the  reference  quan¬ 
tities  with  which  we  form  the  absolute  quantities  used.  Toward 
the  wall  and  indeed  from  ym  =  0.84  (at  A  in  Fig.  2),  we 
then  have  U  =  1.705  -  1.680  ym,  so  that  ym  =  1.015  at  the 
wall  (U  -  0).  Thereby  U  itself  is  represented  with  suffi¬ 
cient  accuracy,  but  we  cannot  put  IP  =  constant  and  Uu  =  0 
in  the.  region  hear  the  wall.  With  ym  as  the  wall  distance 
corresponding  to  the  first  terms  of  the  Blasius  series,  the 
value s 

U’  &  1.68  (1  -  3.65  Yv  )» 

U”  -  -  18.4  y^, 

apply,  instead,  near  the  wall,  ir.  so  far  as  these  quantities 

&P  e  subs  c-ouer.tly 


The  frictionless  solutions  ?,  ana  are  calculated  ac¬ 
cording  to  the  series  formulas  of  Section  II  for  the  parabolic 
distribution  of  U  and  with  the  utilization  of  the  eoluti ohs 


sin  a  y  and  cos  c  y  for  linear  U.  It  was  thu 


aside:  from  the.  irsneiti oral  sursti tut  lc*n,  to  present  the  solu- 
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V 


tj oris  for  the  wall  and  the  parabola  vertex  with  sufficient  accu- 

FT  'i 

racy° ’ -  The  transitional  substitution  cannot,  of  course,  be 

calculated  with  the  roughly  approximated  profile  (straight  line  + 

XJ« 

parabola) ,  but  with  the  more  accurate  value  of  a  -  3,9  c* , 
whibh  is  obtained  from  the  indicated  formulas  for  II B  and  U1 , 
when,  according  to  the  above,  0  m  U  1.68  yw  even  at  the 
critical  poi nt. 

2.  The  intersection  points  of  the  £  curves  (Section  V) 

D 


th  the  -  |p  polar  diagram  are  next  determined  for  0  ~  0. 
An  imaginary  cor:  onent  ie  introduced  into  £  by  the  phasic 

£  is  purely  real  for 


U« 

shift  accompanying  ~  cs,  eo 


m 


sufficiently  small  c.  The  intersection  points  with  the 

diagram  are  therefore  either  at  -  -5-  =  0.58,  or 
.  Tiw 

n 

ox  -  =  0  and  negatively  very  large  fj^. 

only  the  linear  terms  in  c  and  C ,  we  obtain 

-1.1 9  a  +  2o 


1 
% 
=  —  2.3, 


£  i 


0. 598  c  3.36 


Z> 


we 


For  the  first  intersection  point  with  -  Jf.  =  Ov! 

w 

ly  have  o  *  0.74  c  and.,  under  consideration  of  %  =  2.3  or 

(1.68  a  B)a/50.596  c  =  2.3 
R  =  46  c~* 

According  to  equation  (18)  we  have  0  =  1.66  c  for  the  second 

S'A  special  expieesiofc  is  required  by  this  method  for  the  repre¬ 
sentation  of  $;> ,  when  the  critical  point  (U=c)  approaches  very 
near  to  the  junction  point  of  the  straight  line  and  parabola. 
Nevertheless,  re  can  here  dispense  with  the  repetition  of  this 
Consideration,  which  is  important  for  the  numerical  calculation 
but  without  interest  for  the  nob  leu  as  a  whole. 
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-  O-  In  this  region  T*  -Is  very  large 
(2  a  H  c)  s  (1  +  i)  according  to  equation 
i nation  of  E  we  mat  here  include  the 


intersection  with 


higher  terms  in  c  and  a  and  utilize  therefore  the  imaginary 
part  of  D»  to  whioh  the  imaginary  part  of  -  E  should  he 
equal  for  the  transmission  points.  This  becomes  4.36  i  c*  on 


that  is.  fox 


small  c 


tion  of  asymptotic  value  of  D 


between  R,  a  and  c,  both  these  limiting  cases  are  important 
because  bur  asymptotic  calculations  for  large  a  R  apply  to 


them  with  great  accuracy.  The  formulas,  moreover,  yield  a  util 
i sable  approximation  up  to  about  c  =  0.05, 

The  other  transmission  points  are  found  by  the  graphic 
method  developed  in  Section  V  (Fig.  3).  Bevor.d  a  certain  c 
value  Ce.g.  >  in  Fig.  3  at  c  »  0*43),  the  %  curves  yield  ho 
further  ir.  ter  sec  tion  point  with  the  -  .diagram,  though  the; 
otherwise  yield  two  intersection  points.  Before  entering  into 
the  physical  discussion,  we  will  add  a  few  more  remarks  regard 


3.  Our  whole  calculation  is  only  ar  approximation,  which 


s  ooorer  m 


sxitutio 
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0f  Section  lit*  although  the  transitional  srrip  does  not  by  any 
peaiis  always  remain  small*  so  that  we  cart  no  longer  assume 
U|  to  be  approximately  constant  in  the  transitional  strip* 
&nd  even  the  Junction  conditions  have  to  be  changed  ad  regards 
the  greater  jyj  values  on  the  edge  of  the  transitional  strip. 
How,  however,  the  transitional  sub s t i tuti on  can  be  accompli shed 
under  somewhat  different  conditions.  As  the  differential  equa¬ 


tion  for  <f>2 


in  h  in 


,  we  can 


{« fc* 

(19) 


mmsrnm  ii| 


VU*o  22  U*o  SI 


in  which  the  U  derivatives  are  the  original  ones 
y .*)  The  right-hand  member  of  equation  (IS)  is  no 


stant  in  the  transitional  strip. 


can 


con- 

;t 


there  ie  a  solution  of  equation  (IS)  which,  for  large  |h|  , 

Joins  tic  frictionless  solution  <P2  with  the  branch  of  the  log¬ 
arithm  previously  known  as  physically  real.  By  considering  the 
right-hand  i  ember  of  equation  (19)  as  on  inhomogeneous  term  of 
a  differential  equation  of  the  second  order  for  <?”s ,  it  :is 
shown,  with  the  same  means  (asymptotic  discussion  of  Hank  el 
functions  and  of  the  corresponding  integrals)  as  in  Section  III, 
that,  on  the  semicircle  there  introduced,  the  equation 


=  % 


■trn 


tt  j  y 


+  e 


2  a*  (l 


u«0  rn 

4uh2! 


-  -■  at  fr-) 


4;  See  footnote,  page  26. 
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is  identical  with  the  frictionless-disturbanoo  equation,  whereby 
the  transitional  substitution  is  confirmed,  even  under  these 
general  assumptions.  We  commit  :a  slight  error ,  however,  by  tak¬ 
ing  92  at  the  wall  as  already  friotionless  (naturally,  with 
inclusion  of  the  phasic  jump),  i.e. ,  by  assuming  that  the  fric¬ 
tional  effect  has  vanished.  Outside  of  the  transitional  strip 
(for  example,  at  the  point  with  the  index  ta)  the  errors  in  <Pi 

and  <P2,  through  disregard  of  friction  of  the  order  (ct E)~*  , 

• 

are  always  very  small.  On  the  oontrary,  the  error  in  <Pa  ,  of 

0ii 

the  order  e  -—ft,  is  0-055  in  the  most  unfavorable  case, 
b  o 

Tills  rough  estimation  of  the  errors  in  our  approximation 
ehoY/s  that  our  calculation  yields  the  physically  correct  funda¬ 
mental  values,  while  indicating  at  the  same  time  that,  for  the 
last  intersection  points  ( small  Reynolds  Numbers),  the  numerical 
values  of  the  disturbance  parameters  and  of  c  arc  subject  to 
errors  of  a  few  per  cent,  especially  in  the  determination  of 


the  vexoc 


fcxion  o 


or  bare 

25. 

ctional 

terms  <! 

side  of 

equation 

ing  cone 

iceratio 

ffect  tk 

e  t  ran  si 

■ ) 
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smaller  ones,  which  would  have  very  little  effect  on  the  re- 
suits  of  our  calculation. 

4.  In  .giving  the  results  of  our  calculation  of  the  trans¬ 
mission  points,  we  will  write  in  full  the  hitherto  abbreviated 
absolute  quantities.  Instead  of  the  previously  used,  somewhat 
arbitrarily  determined  reference  length  (half  the  width  of  a 
certain  parabola) ,  we  will  introduce  a  physically  logical  ref¬ 
erence  length.  The  determination  of  a  characteristic  length 
from  the  present  and  other  boundary-lay er  profiles  is  rendered 
difficult  by  the  aeymptotio  union  with  the  undisturbed  velocity, 
so  that  nearly  all  of  the  previous  definitions  of  the  boundary- 
layer  thickness  are  not  rational.  The  only  formula  for  the 
boundary-layer  thickness,  which  is  not  arbitrary,  which  can  be 
accurately  determined  by  experiment,  and  which  covers  the  en¬ 
tire  course  of  the  profile  is 

6  =  /  (x  - 

in  which  the  integral  extends  from  the  wall  to  the  undisturbed 
potential  velocity  Um-  This  length,  formerly  termed  also  dis¬ 
placement  thickness  but  here  simply  boundary-layer  thickness* 
yields,  when  multiplied  by  "Ujj.,  just  the  measure  for  the  reduc¬ 
tion  of  the  flow  ty  the  boundary- 1 ay  er  friction.  In  this  case 
6  is  0.341  of  our  former  reference  length.  Vfe  accordingly 
take  as  Reynolds  Number  R. 

According  to  our  determination  of  the  transmission  points, 
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two  disturbance  parameters  correspond  to  each  H,  except  that 
only  one  disturbance  parameter  corresponds  to  the  last  critical 
Reynolds  Humber  430.  Figures  4-5  show  the  disturbance  parame- 


space  enclosed  in  the  curves  is  obviously  the  region  of  the  un- 


that  the  diagrams  cover  only  the  real  part  of  these  quantities 
Here  a  few  transmission  points  are  indicated  numerically 
as  obtained  from  our  calculation  in  connection  with  Figure  3. 
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ey  hardly  continue  to  be  physically  important 
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instability:  ^  ~  0*^25,  0-5  *=  0.367,  ^  =  0.148.  The  great 

*ave  lengths  X  of  the  disturbance  are  especially  striking  in 
coj^-arison  with  the  boundary-layer  thickness*  X  ~  t  sc  that 
jufnirabia  X  =  $  -  17.16.  Also  cex  tain  observations  of 

pr andil  regarding  the  previously  mentioned  channel,  at  the  be¬ 
ginning  of  which  the  flow1  somewhat  re  scales  the  plate  flow,  seem 
to  refer  to  these  great  wave  lengths.  Since  6  increases  along 
the  plate,  it  is  perhaps  desirable  to  assign  to  X  a  limit,  which 


IB  independent  of  x.  According  to  Figure  4,  a  certain  region 

of  unstable  ao  corresponds  to  every  We  can  eliminate  6 

.  tfm6 

from  no  and  and  then  determine  the  maximum  'unstable  a 


or  the  minimum  X  for  any  given 


and  n.  The  minimum  unsta¬ 


ble  X  is  8400  which  is  approximately-proportlonal  to  the 
steady  oscillation  at  R  «  430.  If  the  length  of  the  plate  fror 
its  leading  edge  to  the  point  R  =  =  420  is  designated  by 


1%  =  59000  JLy  because  f  =  1.73  1  /  —iy.  Tr 


X  =  l. 


Tietjens  investigated  the  stability  of  a  profile  which  as¬ 
cends  linearly  to  UK  and  then  continues  with  a  bend.  This 
profile  must  therefore  be  regarded  as  a  rough  approximation  of 
ours.  The  transmission  points  were  taken  from  his  Figure  17 
and  p>lot tod,.,  in  part,  in  o  r  Figure  4  with  crosses,  whereby, 
according  to  Tietjens,  tx  -  2  6 .  Tietjens  obtained  One  branch 
Of  the  equilibrium  curve  in  very  close  approximation.  Since  he 
obtained  .only  true  one  branch,  he  found- no  limit  for'  R  her  for 


'~W\  VJ,*.  J-  remain.  - 
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a.6  which  depends  on  the  general  disappearance  of  the  hend  in 


undertook  to  investigate  the  stability  of  pro 


files  with  bends.  In  general,  Heisenberg  assumes  that  a  ~  0. 
The  convergence  of  the  development  of  the  solutions  according  to 
a  is  not  explained.  The  transitional  substitution  for  q>2 
(Heisenberg  cp4  )  is  accomplished  by  another  and  not  quite  con¬ 
vincing  method.  The  numerical  value  of  the  phasic  bend  a'Tees 
with  ours  in  his  special  case.  Heisenberg  did  not  calculate  a 
profile,  but  confined  himself  to  suppositions  which,  in  part 
(as  regards  the  existence  of  an  upper  limit  fox  a  and  a  lower 
limit  for  R)  lie  in  the  same  direction  as  ovr  results  concern¬ 
ing  rue  special  profile  for  the  flow  along  th<  plate. 

5.  The  comparison  of  our  results  with  those  of  Burgers®) 
and  Eijnen?)  and  Hansen^)  is  incomplete  for  two  reasons.  In  the 
first  place  we  are  not  well  enough  informed  concerning  the  ac¬ 
tual  disturbances,  ana  then  it  is  not  clear  as  to  how  far  the 
turning  point  defined  by  these  writers  agrees  with  our  beginning 
of  instability  of  the  laminar  flow.  Both  Burgers  and  Hansen  ce- 


orwnrd  portion  of 


fine  the  turning  point  as  follows .  Or. 
the  plate  there  is  always  a  laminar  sec 
Frandtl  and  Blasius.  After  the  t 


rang  point,  there  ie  e  decided 
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tuxbulent  portion,  whose  laws  were  explained  by  Prandtl  and  Harman 
(J./7  power  law,  etc*)*  In  both  sections  the  velocity  gradient 
close  to  the  wall  continually  diminishes  downstream  with  an  in- 
crease  only  in  the  transition  zone.  The  critical  Reynolds  Number 
is  calculated  simply  from  the  minimum,  point  of  the  velocity  gra¬ 
dient,  while  the  instability  of  the  laminar  flow  might  begin 
sooner  at  smaller  6  and  hence  at  smaller  R.  At  best  we  may 
therefore  expect  only  a  confirmation  of  our  calculated  critical 
Reynold#  Numbers.  It  should  be  noted  that  the  boundary-layer 
thickness  obtained  by  these  investiga  tors  is  3.18  times  that  ob¬ 
tained  by  us.  On  the  average,  they  obtained  in  our  scale  a  orit- 
ical  Reynolds  Number  of  about  9509^(as  calculated  from  their  val- 
ue  of  3000) ,  which  decreased  with  greater  disturbance  to  50C 
a  disturbing  screen  placed  in  the  air  stream  at  the  leading 
of  the  plate).  If  it  ie  considered  that,  in  our  calculation  at 
420,  there  was  still  a  single  partial  undamped  oscillation,  while , 
in  the  case  of  the  experimentally  determined  turning  point,  there 
was  already  &  considerable  increase  in  the  small  disturbance#,  it 
can  then  be  said,  under  consideration  of  the  above  critical  re¬ 
marks,  that  our  calculated  results  agree  very  well  With  the  ex¬ 
perimental  results, 

r%  \  ~'w-  •~T~'l''ru  '  *  .'-"V"  ■'  * — "  — «— ■'  i  i'  1 1  1 1  ■  1 1  —  ---  -  i  ■  minnm 

Since  our  R  =  ~~  =  1.73  J -  1,73  J where  1  repre- 

6er:x.s  the  plate  length  from  the  leading  edge,  this  numerical 
value  corresponds  to  an  R  •  of  300,000. 
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